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ABSTRACT 

This paper compares recent spatial anomaly time series of OLR (Outgoing Longwave Radiation) and OLR clr (Clear Sky 
OLR) as determined using CERES and AIRS observations over the time period September 2002 through June 2010. We 
find excellent agreement in OLR anomaly time series of both data sets in almost every detail, down to the°lx 1° spatial 
grid point level. This extremely close agreement of OLR anomaly time series derived from observations by two different 
instruments implies that both sets of results must be highly stable. This agreement also validates to some extent the 
anomaly time series of the AIRS derived products used in the computation of the AIRS OLR product. The paper then 
examines anomaly time series of AIRS derived products over the extended time period September 2002 through April 
2011. We show that OLR anomalies during this period are closely in phase with those of an El Nino index, and that 
recent global and tropical mean decreases in OLR and OLR clr are a result of a transition from an El Nino condition at 
the beginning of the data record to La Nina conditions toward the end of the data period. This relationship can be 
explained by temporal changes of the distribution of mid-tropospheric water vapor and cloud cover in two spatial regions 
that are in direct response to El Nino/La Nina activity which occurs outside these spatial regions. 
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1. INTRODUCTION 

OLR (Outgoing Longwave Radiation) is a critical component of the Earth’s radiation budget and represents the total 
radiation going to space emitted by the earth-atmosphere system and integrated over all angles. OLR clr represents the 
Outgoing Longwave Radiation going to space under cloud free conditions. OLR products have been generated and 
monitored globally since 1975 based on broad spectral band measurements taken at a given satellite zenith angle by the 
ERB instrument on the Nimbus-6 and Nimbus-7 satellites, the AVHRR instrument on NOAA operational satellites; and 
most recently by CERES which has flown on EOS Terra since 2000 and EOS Aqua since 2002. 

OLR at a given location is affected primarily by the earth’s skin surface temperature, skin surface spectral emissivity, 
atmospheric vertical temperature profile and water vapor profile, and the heights, amounts, and spectral emissivities of 
multiple layers of cloud cover. OLR can be computed given all the needed geophysical parameters using an OLR 
Radiative Transfer Algorithm (RTA). Mehta and Susskind developed such an OLR RTA used in conjunction with the 
TOVS (TIROS Operational Vertical Sounder) retrieval methodology (Susskind et al. 1997) in order to generate the 
TOVS Pathfinder Path-A OLR data set (Mehta and Susskind 1999a, 1999b). AIRS OLR is computed using AIRS/AMSU 
sounding products in a completely analogous manner, including use of the same Mehta and Susskind OLR RTA 
(Susskind et al. 2003). AIRS OLR clr is also computed for a given scene, otherwise using the same retrieved geophysical 
parameters but setting the amount of clouds to be zero. More details about the AIRS Science Team Version-5 OLR 
products are given in Susskind et al. 2011a and Susskind et al. 2011b. 


2. AIRS AND CERES OLR DATS SETS USED 

In this paper we use the operational monthly mean OLR and OLR clr data products produced by the AIRS and CERES 
Science Teams respectively. We obtained those AIRS OLR products from the Goddard DISC and the CERES products 
from the Langley ASDC. AIRS was launched on the EOS Aqua satellite in a 1:30 AM/PM local crossing time orbit in 



May 2002. The operational processing of AIRS data began after AIRS became stable in September 2002. We use the 
AIRS Version-5 monthly mean Level-3 l°xl° latitude-longitude grid products which contain separate products 
generated for each of the 1:30 AM and PM local time orbits. We averaged the AM and PM products together to generate 
and use a single monthly mean product on the 1° x 1° grid for each month included in the data set. In addition to AIRS 
OLR and OLR clr , we also use the AIRS Level-3 surface skin temperatures, water vapor profiles, and cloud products to 
demonstrate the factors contributing to the anomaly time series of OLR and OLR clr . 

CERES has flown on both EOS Terra, which was launched in December 1999 on a 10:30 AM/PM local crossing time 
orbit, as well as on EOS Aqua, the same platform that carries AIRS. The CERES Science Team generates a number of 
different OLR data sets using CERES observations. The latest versions of the longest record CERES OLR data sets are 
referred to as the CERES SSFldeg-lite Edition-2.5 data sets, which like AIRS, are presented on 1° x 1° latitude-longitude 
grid. CERES SSFldeg-lite Edition-2.5 uses the latest calibration improvements with Edition-2 CERES cloud retrievals 
(Minnis et al. 2008, Minnis et al. 2011), angular dependence models (Loeb et al. 2005), and time-space averaging 
(Young et al. 1998). At the time of this writing, the CERES Terra Edition-2.5 OLR data sets extended to June 2010, the 
CERES Aqua Edition-2.5 OLR data sets extended to February 2010, and the AIRS Level-3 products extended to April 
2011. For this reason, we chose to compare anomaly time series of AIRS OLR and OLR clr with those obtained from 
CERES Terra. 


3. GLOBAL MEAN AND TROPICAL MEAN ANOMALY TIME SERIES 

We generated AIRS and CERES monthly mean OLR and OLR clr climatologies on a 1° x 1° spatial resolution for each 
month by taking the average of the grid box value for that month over a 7-year time period. OLR and OLR clr anomalies 
for a given month in a given year, on a 1° x 1° spatial grid, are defined as the difference between their monthly mean 
values and their monthly climatologies for that grid box. The area mean anomaly for a given month is the cosine latitude 
weighted average of the grid box anomalies contained in the area under consideration. 

Figure la shows the global mean anomaly time series of AIRS Version-5 OLR and CERES Terra OLR for the periods 
September 2002 through April 2011 and September 2002 through June 2010, respectively, as well as the difference 
between the two sets of monthly mean anomalies for the overlap time period. Figure lb shows analogous results for 
tropical mean OLR anomalies, and Figures 2a and 2b show analogous global mean and tropical mean anomaly time 
series for the AIRS and CERES OLR clr products. Figures lb and 2b also include the El Nino index, multiplied by 3 in 
Figure lb and by 2 in Figure 2b. The AIRS El Nino index represents the monthly mean oceanic sea surface temperature 
(SST) anomaly determined using the AIRS Version-5 T skin product, averaged over the spatial area 15°N to 15°S latitude 
and 140°W westward to 160°E longitude. Tropical mean OLR and OLR anomalies tend to track those of the El Nino 
index in phase fairly closely. Amplitudes of the largest tropical mean OLR clr anomalies match those of the El Nino 
index multiplied by 3, and those of OLR clr anomalies match the El Nino index multiplied by 2. This shows that the 
greatest tropical mean OLR anomalies are roughly 50% larger than the greatest tropical mean OLR clr anomalies. 
Positive values of the El Nino index (2003, 2005, 2007, early 2010) correspond to El Nino periods, and negative values 
(2008, late 2010) correspond to La Nina periods. 

Figures la and lb show that an onset of negative global mean OLR anomalies began in late 2007, with tropical mean 
anomalies generally considerably larger than global mean anomalies, especially after mid-2007. The decrease in global 
mean OLR in late 2007 is strongly influenced by the significant reduction in tropical mean OLR which started a few 
months earlier. Tropical mean OLR anomalies, and to a lesser extent global mean OLR anomalies, became positive 
starting in late 2009, roughly coincident with the onset of another El Nino event. The AIRS OLR product, which extends 
to April 2011, shows that very substantial negative global and tropical mean OLR anomalies occurred in the period 
starting mid-2010, in which a substantial La Nina event occurred which continues through the time period observed by 
AIRS shown in this paper. Figures 2a and 2b show analogous features to those in Figure 1. We define the Average Rate 
of Change (ARC) of an anomaly time series as the slope of the linear least squares fit of the time series. The ARC of the 
El Nino index, computed over the extended time period through April 2011, is -O.llK/yr, which indicates that on the 
average, the El Nino region has been cooling over the last 8 years. 
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Fig. la. Global mean AIRS and CERES Terra monthly 
mean OLR anomaly time series and their anomaly 
differences for the overlap period September 2002 through 
June 2010. 
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Fig. lb. Tropical (20°N to 20°S) mean monthly mean OLR 
anomaly time series and anomaly differences of AIRS and 
CERES Terra for the overlap period September 2002 
through June 2010. Also shown in black is the AIRS El 
Nino index multiplied by 3. 
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Fig. 2a. Global mean monthly mean OLR clr anomaly time 
series and anomaly differences of AIRS and CERES Terra 
for the overlap period September 2002 through June 2010. 


Fig. 2b. Tropical (20 °N to 20 °S) mean monthly mean 
OLR clr anomaly time series and anomaly differences of 
AIRS and CERES Terra for the overlap period September 
2002 through June 2010. Also shown in black is the AIRS 
El Nino index multiplied bv 2. 


Table la shows the global and tropical mean values of the ARC’S of AIRS OLR and CERES Terra OLR anomalies over 
the time period September 2002 through June 2010, the standard deviations between the two sets of global and tropical 
anomaly time series, and the temporal correlations between each global and each tropical anomaly time series. The 
agreement of the ARC’S of both global mean and tropical mean anomaly time series found in the AIRS and CERES OLR 
records is very good, with a difference of about 0.02 W/m 2 /yr. The uncertainties shown in this table, and subsequent 
tables, represent the 99% confidence interval. The temporal correlations of the two sets of global and tropical OLR 
anomaly time series are 0.953 and 0.983 respectively. Both show that global mean OLR has decreased on the average on 
the order of-0.075 W/m 2 /yr over the common time period September 2002 through June 2010, and tropical mean OLR 
has decreased at a rate of roughly -0.1 W/m 2 /yr from the beginning of the time period to the end. The close agreement of 
global and tropical mean ARC’S of AIRS and CERES OLR anomaly time series is more significant than the values 
themselves. 

Table lb shows that the AIRS and CERES OLR clr anomaly time series are also highly correlated with each other, 
though the correlation is not as high as between the OLR time series. Part of the decrease in correlation is a result of the 
reduction in values of each set of OLR clr anomalies compared to those of OLR. In addition, the standard deviations of 
the OLR clr anomaly differences are somewhat larger than those of OLR. Perhaps more significant is the fact that AIRS 
and CERES global mean and tropical mean OLR clr anomaly time series are each very highly correlated with those of 
OLR. These and other correlations of anomaly time series are shown in Table lc which shows temporal correlations 
between global mean and tropical mean anomaly time series of OLR and OLR clr as well as correlations of anomaly time 
series with the AIRS El Nino index. Correlations using AIRS time series are shown above the diagonal in bold and those 
using CERES time series are shown beneath the diagonal. The temporal correlation of the CERES global and tropical 
mean OLR anomaly time series is 0.55, and the corresponding correlation for the AIRS time series is 0.56. This shows 
that tropical anomalies provide a significant contribution to the global OLR anomaly time series. The CERES and AIRS 
tropical mean OLR anomaly time series also correlate highly with the El Nino index, with temporal correlations of 0.68 
and 0.69, respectively. The corresponding correlations of the CERES and AIRS OLR clr anomaly time series with the El 
Nino index are 0.70 and 0.63 respectively. Both sets of global OLR and OLR clr anomaly time series also show moderate 
correlations with the El Nino index which are smaller than those found in the tropics. These correlations of global and 
tropical anomaly time series with the El Nino index imply that the recent short term decreases in global and tropical OLR 
are the result of changes from El Nino conditions at the beginning of the time series to La Nina conditions at the end. 


Table la. OLR Anomaly Time Series Comparison 
September 2002 through June 2010 


Data Set 

Global 

Tropical 

AIRS ARC (W/m 2 /yr) 

-0.0881 ±0.0055 

-0.1125 ±0.0128 

CERES Terra ARC (W/m'/yr) 

-0.0655 ± 0.0050 

-0.1006 ±0.0121 

AIRS Minus CERES STD (W/m 2 ) 

0.147 

0.210 

AIRS/CERES Correlation 

0.953 

0.983 


Table lb. OLR clr Anomaly Time Series Comparison 
September 2002 through June 2010 


Data Set 

Global 

Tropical 

AIRS ARC (W/m 2 /yr) 

-0.0045 ± 0.0036 

-0.0177 ±0.0067 

CERES Terra ARC (W/m 2 /yr) 

-0.0703 ± 0.0043 

-0.1018 ±0.0078 

AIRS Minus CERES STD (W/m 2 ) 

0.195 

0.273 

AIRS/CERES Correlation 

0.857 

0.922 





Table lc. Correlations between OLR and OLR clr Anomaly Time Series in Different Domains 

AIRS and CERES 



Global OLR 

Tropical 

OLR 

Global 

OLRclr 

Tropical 

OLRclr 

El Nino 
Index 

Global 

OLR 

— 

0.56 

0.77 

0.48 

0.47 

Tropical 

OLR 

0.55 

— 

0.51 

0.87 

0.69 

Global 

OLR clr 

0.79 

0.56 

— 

0.65 

0.46 

Tropical 

OLR clr 

0.55 

0.87 

0.73 

— 

0.63 

El Nino 
Index 

0.41 

0.68 

0.52 

0.70 

— 


4. THE SPATIAL DISTRIBUTION OF ARC’S 

The spatial distributions of global OLR ARC’S over the time period September 2002 through June 2010 are shown in 
Figures 3a and 3b for AIRS and CERES, respectively. Both fields are presented on the 1° latitude by 1° longitude grid on 
which the data sets are given. Figures 3a and 3b demonstrate two very important points. The first is the virtually 
indistinguishable spatial distributions of the ARC’S of AIRS OLR and of CERES OLR. The spatial correlation of the two 
sets of ARC’S is 0.948, and the standard deviation of their differences is 0.18 W/m 2 /yr. The second is that there is a 
considerable longitudinal variation to OLR ARC’S at a given latitude, particularly in the tropics. While the tropical mean 
ARC of OLR is quite negative, areas of very positive OLR ARC’S exist in the tropics, particularly near the dateline. 

Figures 3a and 3b, and some subsequent figures, contain rectangles surrounding the area between 8°N to 20°S and 140°W 
eastward to 10°E, which we will refer to as OLR Region 1, and the area between 15°S to 30°S and 140°E eastward to 
160°W, which we will refer to as OLR Region 2. ARC’S of OLR within both these areas are very negative over the 
period September 2002 through June 2010. The decreases in OLR that took place in these two regions during the period 
under study will be shown later in the paper to be the major cause of the negative values of the global and tropical mean 
OLR ARC’S shown in Table 1. OLR Regions 1 and 2 were selected entirely based on the results shown in Figure 3 and 
were chosen both so as to be rectangular and also to encompass key features shown in both Figures 3a and 3b. OLR 
Region 2 is essentially in the “heart” of the area referred to as the South Pacific Convergence Zone. 
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Fig. 3a. Spatial 1° latitude by 1° longitude distribution of global 
AIRS OLR ARC’S over the time period September 2002 through 
June 2010. OLR Region 1 is outlined in gray and OLR Region 2 is 
outlined in black in this and most subsequent figures showing 
spatial distributions of ARC’S of different parameters. 


Fig. 3b. Spatial 1° latitude by 1° longitude distribution of 
global CERES OLR ARC’S over the time period 
September 2002 through June 2010 


5. LONGITUDINAL DISTRIBUTION OF EQUATORIAL ANOMALY TIME SERIES: 

HOVMOLLER DIAGRAMS 


Figures 4a and 4b present Hovmoller diagrams, showing time series of monthly mean OLR anomalies (vertical scale), 
integrated over the latitude range 5°N through 5°S, in each 1° longitude bin (horizontal scale) for the periods for which 
AIRS and CERES OLR records exist: September 2002 - April 2011 for AIRS and September 2002 - June 2010 for 
CERES, and Figure 4c shows their difference. The two sets of Hovmoller diagrams show virtually identical features, 
including their magnitudes. The spatial temporal correlation between the two sets of Hovmoller diagrams is 0.993. Most 
of the region covered is ocean. There are three relatively small land areas near the equator: South America, Africa, and 
Indonesia. These land areas each lie between the three sets of narrow vertical lines shown in Figure 4. There are also two 
other gray vertical lines at 140W longitude and 10°E longitude, encompassing the longitudinal range contained within 
Region 1. The anomaly time series shown in Figures 4a and 4b depict the phase relationship of equatorial OLR 
anomalies in different longitudes as a function of time. Such figures in turn provide insight into the spatial distribution of 
tropical ARC’S in the vicinity of the equator shown in Figure 3. In the longitudinal band 161W westward to 140°E, 
equatorial OLR anomalies were very negative in late 2002/early 2003, and were very positive from mid-2007 through 
early 2009. This gave rise to the substantial positive OLR ARC’S shown in Figure 3 over the region 5°N - 5°S, 160°W to 
140°E. Figure 4a shows that AIRS OLR anomalies in this same region are also very positive starting mid-2010. 
Therefore, OLR ARC’S over the extended time period September 2002 - April 2011 would be even more positive in this 
area than those shown in Figure 3a covering a shorter time period. Figure 3 shows very negative values of OLR ARC’S 
near the equator between 100°E and 140°E longitudes. Figure 4 shows that equatorial OLR anomalies between 100°E and 
140°E are out of phase with those between 140°E and 160°W and are of comparable magnitude. Figures 4a and 4b also 
show that equatorial OLR anomalies 160°W eastward to 40°E, within the longitudinal domain of OLR Region 1, tend to 
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Fig. 4a. Hovmoller diagram for time series Fig. 4b. Hovmoller diagram for time Fig. 4c. Hovmoller diagram for the 

of monthly mean AIRS OLR anomalies series of monthly mean CERES OLR difference between AIRS OLR and 

integrated over the latitude range 5°N anomalies integrated over the latitude CERES OLR. 

through 5 °S in each 1° longitude bin for the range 5°N through 5°S in each 1° 

period September 2002 through April longitude bin for the period September 

2011. The longitudinal domain of OLR 2002 through June 2010. 

Region 1 is indicated by the gray vertical 
lines in this and most subsequent 
Hovmoller diagrams. 




be smaller than, and out of phase with, those from 160W westward to 140°E. This gives rise to the negative equatorial 
OLR ARC’S shown in Figure 3 in OLR Region 1. Hovmoller diagrams of AIRS and CERES OLR clr anomalies are very 
similar to each other, with a correlation of 0.913, and are also very similar to those of OLR, but with smaller values of 
the anomalies. 


6. THE EFFECT OF PHASES OF EL NINO/LA NINA ON TROPICAL WATER VAPOR, 
CLOUD COVER, AND OLR ANOMALY TIME SERIES 

Figure 5a shows the spatial distribution of the Average Rate of Change of the AIRS Version 5 surface skin temperature 
T S kin over the period September 2002 - April 2011. A number of important features are found in Figure 5a. While the 
global mean surface skin temperature ARC is essentially zero over this time period, there are areas in which significant 
positive and negative T skin ARC’S exist. Both polar regions have warmed over this time period, especially near the North 
Pole. There is considerable warming and cooling structure in areas over Northern Hemisphere extra-tropical land, and 
there has been substantial cooling over much of Africa, especially south of 15 °S, as well as over much of Australia. The 
tropics are marked by a substantial oceanic cooling over the region 15°N to 15°S and 140°W longitude westward to 
160°E which we refer to as the AIRS El Nino region. The AIRS El Nino region, indicated by the green rectangle in 
Figure 5a, was defined to encompass the region of large negative sea surface temperature ARC’S shown in the figure. 
There has also been substantial oceanic warming over this time period in the areas to the south, west, and north of the 
AIRS El Nino region. The AIRS El Nino region as defined in this paper is similar to the commonly used Nino 4 region, 
5°N - 5°S, 120°W westward to 170°W (Yeh, et al. 2009). 

Figure 5b shows the global distribution of ARC’S of 500 mb specific humidity (%/yr) over the time period September 
2002 to April 2011. OLR is very sensitive to the concentration of mid-upper tropospheric water vapor in very moist (i.e. 
tropical) areas, in the sense that increasing water vapor concentration increases atmospheric absorption in some spectral 
regions and therefore lowers OLR, everything else being equal. The area of strong equatorial surface temperature cooling 
in the AIRS El Nino region between 160°E and 140°W, shown in Figure 5a, is marked by a significant mid-tropospheric 
drying during this period, and the surrounding areas of warming surface skin temperature are marked by considerable 
mid-tropospheric moistening during the period under study. There is also considerable mid-tropospheric moistening from 
8°N to 20° S in the Atlantic Ocean and especially in the Eastern Pacific Ocean over this time period in areas where there 
has been at most a very small increase in surface skin temperature. Figures 5c and 5d show ARC’S of AIRS OLR and 
OLR clr over the extended time period September 2002 - April 2011. Figures 5b-5d all contain the boundaries of OLR 
Regions 1 and 2 defined previously. OLR Regions 1 and 2 both correspond to locations in which there has been 
considerable mid-tropospheric moistening over this time period, as well as considerable decreases in both OLR and 
OLR clr . 

Figures 5c and 5d show that the spatial distribution of ARC’S of AIRS OLR and OLR clr over extra- tropical land areas 
closely track those of T skin with areas of warming skin temperatures corresponding to areas of increasing OLR and 
OLR clr and vice versa for cooling areas. This demonstrates the expected result that a warming surface should result in 
increasing flux to space, all other things being held constant. The spatial relationship between ARC’S of OLR and T skin in 
the tropics is quite different than in the extra-tropics, especially over ocean. There has been a significant equatorial 
cooling between September 2002 and April 2011 in the El Nino region surrounding the equator from 160°E eastward to 
140°W. This is a reflection of the fact that the time period under study started with an El Nino event in 2002, and had La 
Nina events in 2007 and 2010. This change from El Nino conditions to La Nina conditions is the driving force behind the 
recent negative global mean and tropical mean OLR ARC’S. The AIRS El Nino region of T skin cooling is marked by 
significant increases in both OLR and OLR clr . This region is surrounded to the north, south, and west by areas of 
oceanic warming. All of these oceanic areas are marked by large negative OLR ARC’S. It is extremely important to note 
that equatorial OLR also shows large negative ARC’S over OLR Region 1, 8°N to 20°S, 140°W eastward to 10°E, in 
which there are no spatially coherent changes with regard to T skin . In all of these regions, the changes in atmospheric 
water vapor, rather than in surface skin temperature, have been more significant contributors to the ARC’S of AIRS OLR 
and OLR C l R over the extended time period September 2002 - April 2011. 


Figure 5d shows that spatial ARC’S of clear sky OLR, which to first order are independent of changes in cloud cover, 
closely match ARC’S of OLR poleward of°£(both in phase and magnitude. This shows that OLR ARC’S at high 



latitudes are not influenced substantially by changes in cloud cover ae. In the tropics, the patterns of the ARC’S of OLR 
and clear sky OLR are also similar to each other, but the magnitudes of the relative ARC’S are substantially larger with 
regard to OLR than OLR clr . Tropical OLR clr ARC’s are dominated by those of 500 mb specific humidity q 50 o. Tropical 
OLR ARC’s are also affected by changes in specific humidity, but are dominated primarily by those of fractional cloud 
cover ae. ARC’s of fractional cloud cover (not shown) closely match those of 500 mb specific humidity. It is the 
combination of the effects of anomalies of both tropospheric water vapor and fractional cloud cover that give rise to the 
anomalies and ARC’s of OLR as observed by CERES and confirmed by AIRS. 
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Fig. 5a. Spatial distribution of Average Rate of Change of 
AIRS surface skin temperature (K/yr) over the period 
September 2002 through April 2011. The spatial domain of 
the AIRS El Nino region is outlined in green. 


Fig. 5b. Spatial distribution of Average Rate of Change of 
AIRS 500 mb specific humidity (%/yr) over the period 
September 2002 through April 2011. 
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AIRS Version- 5 Clear Sky OLR Anomaly 
Average Rate of Change (Watts/m 2 /yr) 
September 2002 through April 2011 
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Fig. 5c. Spatial distribution of Average Rate of Change of 
AIRS OLR (W/m 2 /yr) over the period September 2002 
through Anril 2011. 


Fig. 5d. Spatial distribution of Average Rate of Change 
of AIRS Clear Sky OLR (W/m 2 /yr) over the period 
September 2002 through April 2011. 



Table 2 shows the spatial correlation of ARC’S of OLR, OLR clr , ae, q 50 o, and T skin over the time period September 2002 
through April 2011. Spatial correlations of ARC’S of these quantities over the latitudinal domain 6tN through 90°N are 
shown above the diagonal in bold, and analogous spatial correlations over the latitudinal domain 15°N through 15°S are 
shown below the diagonal. At high latitudes, shown in bold, ARC’S of OLR and OLR CLR have a spatial correlation of 
0.83 with each other. Both sets of ARC’S are also highly positively correlated with those of T skin , especially OLR clr . At 
high latitudes, ARC’S of specific humidity and fractional cloud cover are not highly correlated with each other or with 
any of the other geophysical parameters shown in the table. 

Table 2. Spatial Correlations of ARC’S of AIRS Derived Products over the Time Period September 2002 through April 

2011 

90°N to 60°N and 15°N to 15°S 



OLR 

OLRclr 

ae 

4500 

T s kin 

OLR 

— 

0.83 

0.07 

0.23 

0.77 

OLRclr 

0.80 

— 

0.43 

0.36 

0.91 

ae 

-0.91 

-0.80 

— 

0.11 

0.44 

4500 

-0.81 

-0.82 

0.69 

— 

0.33 

Tskin 

-0.45 

-0.08 

0.26 

0.43 

— 


Spatial correlations of ARC’S of the same geophysical parameters in the tropics, shown below the diagonal, are very 
different from those at high latitudes. The correlation of ARC’S of OLR and OLR clr in the tropics is again very high 
(0.80). In the North Polar region, the high correlation of ARC’S of OLR and OLR CLR arose because the spatial 
distribution of ARC’S of OLR and OLR clr were both highly positively correlated with those of T skin . In the tropics, this 
high correlation between OLR and OLR clr arises because both are highly negatively correlated with ARC’S of fractional 
cloud cover ae and 500 mb specific humidity q 50 o. Unlike at high latitudes, ARC’S of tropical T skin are not highly 
correlated with those of any other of the geophysical parameters studied. Tropical correlations of ARC’S and OLR and 
OLR clr with those of T skin are actually somewhat negative because tropical correlations of T skin with those of ae and q 50 o 
are positive. The high correlation of ARC’S of OLR clr with those of ae , on which OLR clr does not depend, arises 
because ARC’S of ae and q 50 o are themselves highly positively correlated with each other. 

Figure 6a shows the Hovmoller diagram of monthly mean T skin anomalies for the period September 2002 through April 
2011. The vertical green lines delineate the longitudinal band included in the El Nino region, 160°E longitude eastward 
to 140°W longitude. The largest SST anomalies tend to occur on either side of the dateline between these longitudinal 
limits. Figure 6a demonstrates that the large negative SST ARC near the equator extending from about 160°E eastward to 
140°W is the result of the transition from an El Nino at the end of 2002 to La Ninas over the time periods late 2007 
through 2008, and especially late 2010 through early 2011. Equatorial T skin anomalies between 100°E and 140°E tend to 
be smaller than, and of opposite sign to, those in the vicinity of the dateline. This gives rise to the band of weaker 
positive SST ARC’S near the equator from 100°E to 140°E. 

Figure 6b shows the Hovmoller diagram of 500 mb specific humidity q 50 o. The vertical gray lines delineate the 
longitudinal extent of OLR Region 1, extending from 140°W eastward to 10°E. 500 mb specific humidity anomalies near 
the dateline in general follow those of SST very closely both in magnitude and in phase. This results from the fact that 
positive SST anomalies in the El Nino region correspond to periods of increased convection in that area, leading to 
enhanced moisture of the mid-troposphere. Conversely, negative SST anomalies in the El Nino area correspond to 
periods of decreased convection (increased subsidence) leading to periods of a drier mid- troposphere. Water vapor 
anomalies over Indonesia, from roughly 100°E to 140°E, are not only out of phase with those near the dateline, as are 
those of T skin , but are also of comparable magnitude to water vapor anomalies near the dateline. This is the result of the 
westward shift of the area of maximum convection during La Nina periods from the dateline to Indonesia. This out of 
phase relationship gives rise to the very substantial positive ARC of 500 mb specific humidity over Indonesia during this 
time period, as shown in Figure 5b. Figure 5b also shows substantial positive q 50 o ARC’S in some tropical locations in 




which no significant changes in T skin exist. The most notable of these is off the west coast of South America, in the 
vicinity of 8°N to 20°S from 140°W eastward to 80°W, which is a part of OLR Region 1. There is also another region of 
positive q 50 o ARC’S near the equator going across South America and extending eastward along the Atlantic Ocean to 
about 10°E longitude, also contained within OLR Region 1. Figure 6b shows that equatorial water vapor anomalies off 
the west coast of South America tend to be out of phase with those at the dateline, especially during the large La Nina 
events in 2007-2008 and 2010-2011. This indicates that La Nina periods of decreased convection near the dateline 
correspond to periods of increased convection eastward of 140E. The same relationship is found to a lesser extent over 
the Atlantic Ocean extending to 1 0°E longitude at the eastern end of OLR Region 1 . 


Figure 6c shows the Hovmoller diagram of the AIRS effective cloud fraction ae. The anomalies of effective cloud 
fraction averaged over °N to 5°S latitude are closely in phase with those of 500 mb specific humidity, with a 
temporal/spatial correlation of 0.717. This correlation is a result of the fact that periods of increased convection (or 
subsidence) correspond to increases (or decreases) of not only mid-tropospheric water vapor but also cloud cover, 
especially with regard to high clouds. Anomalies of OLR clr are driven primarily by those of q 50 o, and those of OLR are 
driven primarily by effective cloud fraction, ae , with a smaller contribution from q 50 o. The high correlation between 
tropical anomalies of OLR and OLR clr are the result of the correlation between anomalies of ae and q 50 o. 
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Fig. 6a. Hovmoller diagram of AIRS 
Surface Skin Temperature anomalies 
(K). The longitudinal domain of the El 
Nino region is enclosed between the 
green vertical lines. 


Fig. 6b. Hovmoller diagram of AIRS 500 
mb specific humidity anomalies (%). The 
longitudinal domain of Region 1 is 

enclosed between the gray vertical lines. 


Fig. 6c. Hovmoller diagram of AIRS 
effective cloud fraction anomalies 
(%). The longitudinal domain of 
Region 1 is enclosed between the 
gray vertical lines. 


Table 3 shows the spatial-temporal correlations of the Hovmoller diagrams of anomalies of OLR, OLR clr , ae, q 50 o, and 
T skin derived from AIRS observations over the period September 2002 through April 2011. Unlike Table 2, this table is 
symmetric above and below the diagonal because only one set of correlation statistics is presented. Spatial-temporal 
anomaly correlations between OLR, OLR clr , ae, q 50 o, and T skin averaged over the latitude range !N through 5°S are 
extremely close to those shown in the lower half of Table 2, representing spatial correlations of ARC’S of these 
geophysical parameters in the tropics. This shows that not only are the spatial patterns of Average Rates of Change of 
tropical OLR, OLR clr , ae, and q 50 o highly correlated with each other, but the detailed structure of their spatial anomalies 
as a function of time are highly correlated as well. 


Table 3. Spatial-temporal Correlations of AIRS Anomaly Hovmoller Diagrams 5N -5°S September 2002 
Through April 2011 
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7. ATTRIBUTION OF RECENT DECREASES IN GLOBAL AND TROPICAL MEAN OLR 
TO CHANGES CONTAINED WITHIN OLR REGIONS 1 AND 2 


Figure 7a shows the AIRS OLR anomaly time series averaged over OLR Region 1 in dark blue. Figure 7a also includes 
in black the El Nino index, this time multiplied by 4, rather than multiplied by 3 as was done in Figure lb. The Region 1 
OLR anomaly time series closely follows that of the El Nino index, but is lagged in time (El Nino occurs first) by about 3 
months. Therefore, in Figure 7a, Region 1 OLR anomalies are plotted 3 months earlier than they actually occurred. For 
this reason, the dark blue line starts 3 months before the black line, and also ends 3 months earlier than the black line. 
The correlation coefficient of the 3 month lagged Region 1 OLR anomaly with that of the El Nino index is 0.898. Figure 
7a also contains the anomaly time series of OLR clr averaged over Region 1, shown by the light blue line, which like the 
dark blue line is lagged 3 months from the El Nino index. OLR clr anomalies in Region 1 are smaller than those of OLR, 
but are highly correlated with those of OLR, with a correlation coefficient of 0.873. The 3 month lagged OLR clr 
anomalies are also highly correlated with those of the El Nino index, with a correlation of 0.821. 


AIRS Version-5 Regional Anomaly Time Series 


September 2002 through April 2011 



OLR averaged over Region 1 lagged by 3 months 


OLR clr averaged over Region 1 lagged by 3 months 

Fig. 7a. AIRS OLR and OLR clr anomaly time series averaged 
over OLR Region 1 superimposed on the AIRS El Nino index 

multiplied by 4 and plotted 3 months earlier. 


AIRS Version-5 Regional Anomaly Time Series 
September 2002 through April 2011 



OLR averaged over Region 2 
OLR clr averaged over Region 2 

Fig. 7b. AIRS OLR and OLR clr anomaly time series 
averaged over OLR Region 2 superimposed on the 
AIRS El Nino index multiplied by 10 without a time lag. 


The OLR and OLR clr anomalies averaged over OLR Region 1 are considerably larger than corresponding tropical mean 
anomalies shown in Figures lb and 2b. Consequently the negative ARC’S of OLR and OLR clr averaged over OLR 
Region 1, which are -0.530 W/m 2 yr and -0.150 W/m 2 /yr respectively, are considerably larger than the corresponding 
negative ARC’S of tropical mean OLR and OLR clr given by -0.188 W/m 2 /yr and -0.074 W/m 2 /yr. These values are 
included in Table 4, which shows area mean Average Rates of Change of the AIRS anomaly time series of OLR and 
OLR clr , computed over the time period September 2002 through April 2011, averaged over different spatial domains. 
The negative tropical mean OLR ARC over the period September 2002 through April 2011 computed as previously, but 
after replacing OLR ARC’S in OLR Region 1 by zeroes, is reduced from -0.188 W/m 2 /yr to -0.032 W/m 2 /yr, and the 
negative global mean OLR ARC is reduced from -0.093 W/m 2 /yr to -0.039 W/m 2 /yr. When an analogous procedure is 
done with regard to OLR clr , the negative tropical mean OLR clr ARC is reduced from -0.074 W/m 2 /yr to -0.030 
W/m 2 /yr and the negative global mean ARC of OLR clr is reduced from -0.018 W/m 2 /yr to -0.003 W/m 2 /yr. This 
indicates that a large part of the recent negative global mean and tropical mean OLR and OLR clr ARC’s results from the 
contribution of OLR and OLR clr anomalies within OLR Region 1. These anomalies are in turn in phase with the El Nino 
index, but lagged in time by roughly 3 months. It is for this reason that the global and tropical OLR and OLR clr anomaly 
time-series shown in Figures 1 and Figure 2 are in phase with, and highly correlated with, El Nino/La Nina activity as 
noted earlier. 

Figures 5c and 5d also show that there are substantial negative ARC’s of OLR and OLR clr within the rectangular box 
surrounding the area from 15°S to 30°S and westward from 150°W to 140°E, referred to as OLR Region 2. One-third of 
this region is in the tropics and the remaining two-thirds is in the subtropics. OLR Region 2 is also marked by large 
positive ARC’s of mid-tropospheric water vapor, as shown in Figure 5b, and cloud cover (not shown) over the period 
under study. Figure 7b shows in dark blue and light blue the time series of OLR anomalies and OLR clr anomalies in 
OLR Region 2, superimposed on the time series of the El Nino index, this time multiplied by 10, without a time lag. As 
shown in Table 4, OLR in Region 2 has an extremely large negative ARC of -1.74 W/m 2 /yr, and OLR clr in Region 2 has 
a very large negative ARC of -0.54 W/m 2 /yr. All three curves in Figure 7b are highly correlated with each other in time. 
The correlation between the unlagged anomaly time series of the El Nino index and OLR in region 2 is 0.854, the 
correlation between the El Nino index and the anomaly time series of OLR clr in Region 2 is 0.781, and the correlation 
between the anomaly time series of OLR and OLR CLR in region 2 is 0.936. As shown in Table 4, the ARC’s of global and 
tropical mean OLR and OLR clr , computed after zeroing out OLR and OLR clr ARC’s in both OLR Region 1 and OLR 
Region 2, are reduced to essentially zero. This shows that the large negative global and tropical ARC’s of OLR and 
OLR clr over the period September 2002 through April 2011 can be attributed almost completely to temporal changes in 
mid-tropospheric water vapor and cloud cover that have taken place primarily within OLR Region 1 along the equator in 
the eastern Pacific and Atlantic Ocean, with a further contribution from OLR Region 2 primarily in the southern 
hemisphere subtropics south-southwest of the El Nino region, covering Eastern Australia and parts of the southwestern 
Pacific Ocean. 


Table 4. Area Mean Average Rates of Change of OLR and OLR clr (W/m 2 /yr) 
September 2002 through April 2011 


Spatial Area 

OLR 

OLRclr 

Global 

-0.093 ± 0.0045 

-0.018 ±0.0031 

Tropical 

-0.188 ± 0.0111 

-0.074 ± 0.0059 

Region 1 

-0.530 ± 0.0229 

-0.150 ±0.0097 

Region 2 

-1.740 ±0.0672 

-0.540 ±0.0218 

Global excluding Region 1 

-0.039 ± 0.0035 

-0.003 ± 0.0026 

Tropical excluding Region 1 

-0.032 ± 0.0080 

-0.030 ±0.0038 

Global excluding Region 1 and 2 

-0.004 ± 0.0035 

0.008 ± 0.0026 

Tropical excluding Region 1 and 2 

0.002 ± 0.0078 

-0.018 ±0.0035 
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